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We report first results on the calculation of NNLO corrections to event shape distributions in 
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corrections yields a considerably better agreement between theory and experimental data both in 
shape and normalisation of the event shape distributions. 



8th International Symposium on Radiative Corrections (RADCOR) 
October 1-5 2007 
Florence, Italy 



* Speaker. 

^Supported in part by the Swiss National Science Foundation (SNF) under contract 200020-117602. 
■'"Supported by the UK Science and Technology Facilities Council. 



© Copyright owned by the author(s) under the terms of the Creative Commons Attribution-NonCommercial-ShareAlike Licence. 



http://pos.sissa.it/ 



e+e event shapes at NNLO T. Gehrmann 



1. Introduction 

For more than a decade experiments at LEP (CERN) and SLC (SLAC) gathered a wealth of 
high precision high energy hadronic data from electron-positron annihilation at a range of centre- 
of-mass energies [[I], |2|, ||]. This data provides one of the cleanest ways of probing our quantitative 
understanding of QCD. This is particularly so because the strong interactions occur only in the fi- 
nal state and are not entangled with the parton density functions associated with beams of hadrons. 
As the understanding of the strong interaction, and the capability of making more precise theo- 
retical predictions, develops, more and more stringent comparisons of theory and experiment are 
possible, leading to improved measurements of fundamental quantities such as the strong coupling 
constant [J]. 

In addition to measuring multi-jet production rates, more specific information about the topol- 
ogy of the events can be extracted. To this end, many variables have been introduced which charac- 
terise the hadronic structure of an event. With the precision data from LEP and SLC, experimental 
distributions for such event shape variables have been extensively studied and have been compared 
with theoretical calculations based on next-to-leading order (NLO) parton-level event generator 
programs [^, |6|, 0|, improved by resumming kinematically-dominant leading and next-to-leading 
logarithms (NLO-i-NLL) ^ and by the inclusion of non-perturbative models of power-suppressed 
hadronisation effects [^. 

The precision of the strong coupling constant determined from event shape data has been lim- 
ited up to now largely by the scale uncertainty of the perturbative NLO calculation. We report 
here on the first calculation of NNLO corrections to event shape variables, and discuss their phe- 
nomenological impact. 

2. Event shape variables 

In order to characterise hadronic final states in electron-positron annihilation, a variety of event 
shape variables have been proposed in the literature, for a review see e.g. [p^]. These variables 
can be categorised into different classes, according to the minimal number of final-state particles 
required for them to be non-vanishing: In the following we shall only consider three particle final 
states which are thus closely related to three-jet final states. 



Among those shape variables, six variables [11| were studied in great detail: the thrust T, the 



normalised heavy jet mass p, the wide and total jet broadenings and Bj, the C-parameter and 
the transition from three-jet to two-jet final states in the Durham jet algorithm F3. 

The perturbative expansion for the distribution of a generic observable y up to NNLO at e+e^ 
centre-of-mass energy y/s, for a renormalisation scale }j} is given by 



1 = ffM*Q)^d+fMi!)Vf^+^ftiog''' 



Chad \ 2k J dy \ 2k J \dy Ay s 



+^(j82log2^+^ilog^^) ]+&{aA). (2.1) 
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The dimensionless perturbative coefficients A, B and C depend only on the event shape variable y. 
They are computed by a fixed-order parton-level calculation, which includes final states with three 
partons at LO, up to four partons at NLO and up to five partons at NNLO. LO and NLO corrections 
to event shapes have been available already for a long time [§, ^, 0]. 

The calculation of the NNLO corrections is carried out using a newly developed parton-level 
event generator programme EE RAD 3 which contains the relevant matrix elements with up to five 



external partons [|12[ [13|, |14[ [15|]. Besides explicit infrared divergences from the loop integrals, 
the four-parton and five-parton contributions yield infrared divergent contributions if one or two of 
the final state partons become collinear or soft. In order to extract these infrared divergences and 
combine them with the virtual corrections, the antenna subtraction method [ p^ was extended to 



NNLO level [17] and implemented for e^e^ 3 jets and related event-shape variables [18]. The 
analytical cancellation of all infrared divergences serves as a very strong check on the implemen- 
tation. EERAD3 yields the perturbative A, B and C coefficients as histograms for all infrared-safe 
event-shape variables related to three-particle final states at leading order. From these. A, S and 
C are computed by normalising to the total hadronic cross section. As a cross check, the A and 
B coefficients have also been obtained from an independent integration [Q] of the NLO matrix 
elements showing excellent agreement. 



For small values of the event shape variable y, the fixed-order expansion, eq. (gjj), fails to 
converge, because the fixed-order coefficients are enhanced by powers of ln(l/y). In order to 
obtain reliable predictions in the region of y ^ 1 it is necessary to resum entire sets of logarithmic 
terms at all orders in a^. A detailed description of the predictions at next-to-leading-logarithmic 



approximation (NLLA) can be found in Ref. [19]. 
3. NNLO results 

The precise size and shape of the NNLO corrections depend on the observable in question. 
Common to all observables is the divergent behaviour of the fixed-order prediction in the two-jet 
limit, where soft-gluon effects at all orders become important, and where resummation is needed. 
For several event shape variables (especially T and C) the full kinematical range is not yet realised 
for three partons, but attained only in the multi-jet limit. In this case, the fixed-order description 
is also insufficient since it is limited to a fixed multiplicity (five partons at NNLO). Consequently, 
the fixed-order description is expected to be reliable in a restricted interval bounded by the two-jet 
limit on one side and the multi-jet limit on the other side. 

In this intermediate region, we observe that inclusion of NNLO corrections (evaluated at the Z- 
boson mass, and for fixed value of the strong coupling constant) typically increases the previously 
available NLO prediction. The magnitude of this increase differs considerably between different 
observables[||], it is substantial for T (18%), Bj (17%) and C (15%), moderate for p and Bw 
(both 10%) and small for F3 (6%). For all shape variables, we observe that the renormalisation 
scale uncertainty of the NNLO prediction is reduced by a factor 2 or more compared to the NLO 
prediction. Inclusion of the NNLO corrections modifies the shape of the event shape distributions. 
We observe that the NNLO prediction describes the shape of the measured event shape distributions 
over a wider kinematical range than the NLO prediction, both towards the two-jet and the multi-jet 
limit. To illustrate the impact of the NNLO corrections, we compare the fixed-order predictions for 
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Figure 1: Perturbative fixed-order predictions for the is -distribution, compared to LEP2 data from ALEPH. 



to LEP2-data obtained by the ALPEH experiment in Figure [IJ, which illustrates especially the 
improvement in the approach to the two-jet region (large — ln(F3)). 

The information contained in the event shape distributions can be restructured by comput- 
ing individual moments. Moments of event shape distributions have been studied theoretically 
and experimentally in particular in view of understanding non-perturbative power corrections |^]. 
Consequently, perturbative NNLO corrections will improve the discrimination between higher per- 
turbative orders and genuine non-perturbative effects. For the first moment {I — T) of the thrust 
distribution, we find the integrated coefficients 

= 2.101 =:^ = 44.98 ^= 1095 ±130, 

which yields for^/s = |J.= Mz'. 

(l-r)(a,(Mz) =0.1189) =0.0398 (LO) + 0.0146(NLO) + 0.0068 (NNLO) . 

Work on moments of the event shapes is ongoing. 
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4. Determination of the strong coupling constant 

Using the newly computed NNLO corrections to event shape variables, we performed[^T]] a 
new extraction of a, from data on the standard set of six event shape variables, measured by the 
ALEPH collaboration [|l|] at centre-of-mass energies of 91.2, 133, 161, 172, 183, 189, 200 and 206 
GeV. The combination of all NNLO determinations from all shape variables yields 

a, (Mz) =0.1240 ± 0.0008 (stat) ± 0.0010 (exp) ± 0.0011 (had) ± 0.0029 (theo) . 

We observe a clear improvement in the fit quality when going to NNLO accuracy. Compared to 
NLO the value of is lowered by about 10%, but still higher than for NLO+NLLA [p, which 
shows the obvious need for a matching of NNLO+NLLA for a fully reliable result. The scatter 
among the ttv-values extracted from different shape variables is lowered considerably, and the 
theoretical uncertainty is decreased by a factor 2 (1.3) compared to NLO (NLO+NNLA). 

These observations visibly illustrate the improvements gained from the inclusion of the NNLO 
corrections, and highlight the need for further studies on the matching of NNLO+NLLA, and on 
the derivation of NNLLA resummation terms. 

5. Outlook 

Our results for the NNLO corrections open up a whole new range of possible comparisons 
with the LEP data. The potential of these studies is illustrated by the new determination of re- 
ported here, which can be further improved by the matching NLLA+NNLO, currently in progress. 
Similarly, our results will also allow a renewed study of power corrections, now matched to NNLO. 
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